function ejer6a
    
    A = [3 1 1;
		 2  6 1;
		 1 1  4] 
    
    b = [5;
		 9;
		 6]
    
    x = [0;
		 0;
		 0];
		 
	er = 10^-11;
	m = 1000;
	
	fprintf('solucion de octave \n');
    A\b
    
    fprintf('solucion de jacobi \n');
    [x1,k1] = jacobi(A, b, er, m, x)
    
    fprintf('solucion de gauss seidel \n');
    [x2,k2] = gseidel(A, b, er, m, x)
    
    
    
 endfunction   